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PRODUCT CONSTRUCTION OF AFFINE CODES* 

YEOW MENG CHEE+, HAN MAO KIAH+, PUNARBASU PURKAYASTHAt, AND 
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Abstract. Binary matrix codes with restricted row and column weights are a desirable method 
of coded modulation for power line communication. In this work, we construct such matrix codes that 
are obtained as products of affine codes - cosets of binary linear codes. Additionally, the constructions 
have the property that they are systematic. Subsequently, we generalize our construction to irregular 
product of affine codes, where the component codes are affine codes of different rates. 
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1. Introduction. Product codes were introduced by Elias [10] and subsequently 
generalized by Forney [11] to concatenated codes. Product codes are a method of 
constructing larger codes from smaller codes while retaining the good rates and good 
decoding complexity from the smaller codes. The codewords of a product code can 
be written as matrices with the rows belonging to the row component code and the 
columns belonging to the column component code. List decoding algorithms have 
also been studied in this context in Barg and Zemor [4] where the min-sum algorithm 
was shown to be amenable to list decoding of product codes. 

Product codes have been subsequently generalized to yield codes obtained from 
product of nonlinear codes by Amrani [2], and to multilevel product codes by Zinoviev 
[19]. Amrani [2] gave the construction of product codes from component nonlinear 
codes which are binary and systematic. The construction guarantees that all the 
columns of any codeword belong to the column component code; however only the 
first few rows corresponding to the systematic part of the column code are guaranteed 
to belong to the row code. In the case where one of the component codes is linear, 
Amrani [2] proposed two soft-decision decoding algorithms. Irregular product codes, 
introduced by Alipour et al. [1], are yet another generalization of product codes where 
each row and column code can be a code of different rate. Irregular product codes 
were introduced to address the need for unequal error protection from bursty noise 
when some parts of the codeword are more vulnerable to burst errors than others. 

In this work we study constructions of systematic nonlinear product codes which 
are obtained as products of affine codes - cosets of linear codes. In contrast to the work 
of Amrani [2], our construction guarantees that all the rows belong to the (affine) row 
code and all the columns belong to the (affine) column code. One primary motivation 
for studying such class of codes arises from a previous study on coded modulation 
for power line channels by Chee et al. [7] that proposed a generalization of the coded 
modulation scheme of Vinck [18]. 

Chee et al. [7] showed that binary matrix codes with bounded column weights, in 


*This paper was presented in part at IEEE International Symposium on Information Theory 2014. 
Submitted: 10 July 2014. Accepted: 03 June 2015. 

t Division of Mathematical Sciences, School of Physical and Mathematical Sci¬ 
ences, Nanyang Technological University, Singapore (ymchee _at_ ntu _dot_ edu _dot_ sg, 
hmkiah _at_ ntu _dot_ edu _dot_ sg) 

*CGG Services Pte Ltd, Singapore (punarbasu+paper _at_ gmail _dot_ com) 

§Institut Mines-Telecom, Telecom ParisTech, CNRS LTCI, France 
(patrick.sole _at_ telecom-paristech _dot_ fr) 


1 



2 


Y. M. CHEE, H. M. KIAH, P. PURKAYASTHA, P. SOLE 


conjunction with multitone frequency shift keying, can be used to counter the harsh 
noise characteristics of the power line channel. Concatenated codes obtained from 
the concatenation of constant weight inner codes with Reed-Solomon outer codes 
were used to obtain families of efficiently decodable codes with good rates and good 
relative distances. In this work, we continue this line of investigation and introduce 
binary systematic product codes with the additional restriction that the row and 
column weights are bounded. The restriction on the column weight arises from the 
desire to be able to detect and correct impulse noise that is present in the power 
line channel. The restriction on the row weights allows one to detect and correct 
narrowband noise. It is quite evident that product codes obtained from the product 
of linear codes do not satisfy these restrictions. The nonlinear codes studied in this 
paper are constructed to satisfy these properties. The efficient decoding algorithms of 
product codes are directly applicable to the constructions presented in this work. As 
a first step to the decoding process, we subtract the coset representative that is used 
in the construction. The coset representative is explicitly described, as explained in 
the following sections. 

The rest of the paper is organized as follows. In the next section we introduce the 
basic definitions and notation that are used throughout the rest of this paper. Section 
3 discusses the general construction of q- ary systematic codes which are product of 
affine codes. Section 4 uses the construction in Section 3 to give constructions of binary 
product codes with restricted row and column weights. This is of interest because of 
its application to coded modulation for power line channels. In Section 5, we extend 
this construction to product codes which can provide unequal error protection, where 
different rows and columns belong to different row and column codes. This section 
generalizes the irregular product code construction of Alipour et al. [1], where the 
component codes are linear codes, to irregular product codes where the component 
codes are affine codes. 

2. Notation and Definitions. Denote the set of integers {1,2,...,, n} by [n] 
for a positive integer n. Denote the finite field of order q by F 9 . A code C of length 
n is a subset of F”, while a linear code C of length n is a linear subspace of F”. The 
dimension of a linear code C is given by the dimension of C as a linear subspace of 
F{h Elements of C are called codewords. Endow the space F” with the Hamming 
distance metric and for u£FJ, the Hamming weight of u is the distance of u from 
the all-zero codeword. A code C C F£ is said to have distance d if the (Hamming) 
distance between any two distinct codewords of C is at least d. Moreover, a linear 
code C has distance d if the weight of all nonzero codewords in C is at least d. We use 
the notation [n, k , d] to denote a linear code of length n, dimension k and distance d. 

Let m, n be positive integers and let F™ xrl denote the set of m by n matrices 
over F g . The transpose of a matrix M is denoted by M T and we regard the vector 
u E F™ as a row vector, or a matrix u in Fj x ". Hence, u 1 denotes a column vector 
in F™ xl . In addition, let 0„ and j n denote the all-zero and all-one vectors of length 
n respectively, while I n and 0 mxn denote the ( n x n ) identity and (to x n ) all-zero 
matrix respectively. We denote the span of a vector u by the notation (u). 

Let C be a linear [ n , k, d\ code. After a permutation of coordinates, there exists a 
matrix A E ¥q X ^ n k ^ such that each codeword in C can be written as (x, xA), where 
x E Fq is called the information vector. The matrix (I/-|A) is said to be a systematic 
encoder of C. 

Let Ci and Ci be linear [n,ki,di] and [n, A^cfe] codes, respectively. Suppose 
Ci C C 2 and pick ueC 2 . Then the set of codewords Ci + u is a coset of Ci in C 2 ■ The 
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collection of all cosets of C\ in C 2 is denoted by C 2 /Ci. Moreover, any coset in C 2 /Ci 
is a (n,di) code of size q kl , and we call the coset an affine [n, ki,di) code. 

Observe that if (1^ |Ai) is a systematic encoder for C\ and u = (ui, u 2 ) where Ui 
is of length ki, then C\ + u = Ci + (Ofe 15 u 2 — uiAi). On the other hand, every coset 
in C 2 /C 1 contains at most one element of the form (O^a). Hence, for every coset 
Ci + u, there is exactly one element of the form (CL,, a), and in this paper, we refer to 
this element as the coset representative of Ci + u. The set of all coset representatives 
of cosets in C 2 /Ci is denoted (C 2 /Ci) rep . 

We also consider the notion of systematicity for nonlinear codes. Let C be a 
(matrix) code of size q k . Then C is said to be systematic of dimension k if there exists 
k coordinates such that C when restricted to these k coordinates is F£. Observe that 
if C is a linear [n, fc, d] code, then any affine code C + u is systematic of dimension k. 

2.1. Matrix Codes. An (to x n)-matrix code C is a subset of F™ xn , while a 
linear (to x n)-matrix code C is a linear subspace of F™ x ", when considered as a vector 
space of dimension mn. Regarding each matrix in F™ x " as a vector of length mn, we 
have the definitions of Hamming distance, Hamming weight and dimension. A linear 
(to x n)-matrix code of dimension K and distance d is denoted by [to x n, K, d\. 

2.2. Classical Product Codes. The classical product code constructs matrix 
codes from two linear codes. Given a linear [n,k,dc] code C and a linear [to, Z,dx>] 
code V, let (Ifc|A) and (1/ |B) be their respective systematic encoders. The product 
code , denoted by C <S> D, is then given by the (to x n)-matrix code (see [14, p. 568]) 


C<8>£> = 


( 


M 

MA 

^ B t M 

b t ma 


: M e F^ xfc , 


where M corresponds to the information bits. It can be shown that C ® T> is a linear 
[to x n, kl, dx>dc\ code. Furthermore, C <g> V has the following property that depends 
on the component codes C and T>. 

Property (C, V). For every N g C (g> V, 

(i) every row of N belongs to C, and 

(ii) every column of N belongs to T>. 

In this paper, we consider nonlinear component codes. Specifically, let C be a 
nonlinear code of length n and size q k and D' be a nonlinear code of length to and 
size q l . We aim to construct an (m x n)-matrix code C ® V of size q kl such that 
Property (C',V r ) holds. This construction differs from the nonlinear product code 
construction in Arnrani [2] because we guarantee that all the rows in every codeword 
belong to the row code C. 

3. Product Codes from Affine Codes. In this section, we provide the general 
construction of systematic matrix codes that are obtained as products of cosets of 
linear codes, i.e., as products of affine codes. Throughout this section, let C and V be 
a linear [n,k,dc] an d [to, l,dx>\ codes respectively. We consider affine codes that are 
obtained as cosets of the codes C and V, i.e., they are of the form C + u and V + v, 
respectively, where u and v are of lengths n and to respectively. In particular, we show 
that if both C and T> contain the all-one vector, then there exists an (to x n)-matrix 
code, that is systematic of dimension kl with Property (C + u, D + v). 
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(3.1) 

(C + u)®(P + v) 4 

(3.2) 

(C + u)®(2? + v) 4 


/ 

V 


M 

MA +j[sl 

B T M + b T j fc 

(B t M + b T j fc )A + jj n _ l a 


: M G F' 


Ixk 


M 

MA+jfa ^ 

v B t M + b T j fe 

B T (MA + jfa) + b T j n _ fc ^ 


: M G F' 


Ixk 


Let (Ifc|A) and (1/ |B) be systematic encoders for C and V respectively. Recall 
that the set of coset representatives of cosets of C\ in C 2 is denoted by (C 2 /Ci) rep . 
Without loss of generality, pick u = (Ofc,a) G (F™/C) rep and v = (0;,b) G (F™/2?) rep . 
Then a typical element in C + u is of the form (x, xA + a) where x is the information 
vector of length k. Similarly, a typical element in V + v is of the form (x, xB + b) 
where x is the information vector of length l. 

Define (C + u) ® (T> + v) to be the (to x n)-matrix code given by (3.1). This is 
obtained by the encoding the first k columns by V + v, followed by encoding all the 
rows by C + u. We observe that for every N G (C + u) ® (T> + v), each row of N 
belongs to C + u. However, we can guarantee only that the first k columns belong to 
V + v. 

On the other hand, if we alter the definition given in (3.1) to be (3.2), where we 
encode the first l rows by C + u, followed by encoding all the columns using V + v, 
we have that every column of N belongs to V + v for each N G (C + u)(g>(£> + v). 

Therefore, the matrix code (C + u) (g> (T> + v) meets our requirements if 

(B t M + b T jfc)A + jm-i a = B t (MA +j? a) + b T j n _ fc ,that is, 

(3.3) b r (j fc A - j n _fc) = (B T jf - j„_i) a 


If (3.3) holds, then (C + u) ® (£> + v) (or equivalently, (C + u)®(Z? + v)) is a coset 
of C <g> D. That is, (C + u)®(D + v) = (C®P) + U, where 


(3.4) 


A 

O 

X 

?r 

jla ' 


v b T j fe 

b T j fe A+j^_ ; a t 


Theorem 3.1. Let C and D be linear [n,k,dc\ and [to, l,dx>\ codes respectively 
and (Ifc|A) and (I;|B) be their respective systematic encoders. Pick u = (Ojt,a) G 
(F£/C) rep and v = (0,,b) G (F™/2?) rep . 

If in addition (3.3) holds, then (C + u) ® (T> + v) defined by (3.1) is equal to (C + 
u)<8>(r> + v) defined by (3.2). Moreover, the code is systematic of dimension kl and 
is a coset of C ®T> with Property (C + u, 2? + v). 

We now provide a sufficient condition for (3.3) to hold. Observe that j„ G C 
if and only if j^.A = j since jfc(Ifc[A) is necessarily j n . This is because the 
only message vector that can give rise to the all-one vector must have all-one in the 
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systematic part of the codeword. Hence, j^.A — j n -k = On-k and b r (jj,A — j„_fc) = 
0 (m-i)x(n-k)- Similar argument holds for B T j J — j ^ n _ l . Hence, (3.3) holds and the 
coset representative U given by (3.4) is 


(3.5) 


( 

T \ 

0/x k 

j/ a 

K b T j fc 

b T j„- fc +j^_ z a y 


and is independent of the matrices A and B. The following corollary, that we refer 
to as Construction J, is now immediate. 

Corollary 3.2 ( Construction I). Let C and T > be linear [n, k, dc] and [to, Z,dx>] 
codes respectively and (Ifc|A) and (I/|B) be their respective systematic encoders. Pick 
u = (Ofc,a) 6 (F^/C) rep and v = (0;,b) G (F™/2?) rep . If in addition j„ G C and 
j m G T>, then (C + u) ® (D + v) defined by (3.1) is systematic of dimension kl and is 
a coset of C ®T> with Property (C + u, V + v). 

Binary linear codes that contain the all-one vector are called self-complementary 
codes. Well-known examples of linear self-complementary codes include the primitive 
narrow-sense Bose-Chaudhuri-Hocquenghem codes, the extended Golay code and the 
Reed-Muller codes [14] . Examples of g-ary linear codes that contain the all-one vector 
include the Reed-Solomon codes, generalized Reed-Muller codes, [14], and difference 
matrix codes [5]. 

4. Variants of Construction I. In this section, we adopt Construction I to 
certain nonlinear component codes C', V that are variants of cosets of linear codes. 
Several well-known families of nonlinear codes, such as Nordstrom-Robinson, Delsarte- 
Goethals, Kerdock and Preparata, can be obtained as unions of cosets of linear codes 
(see [14, Ch. 15]). In general, it is difficult to achieve a matrix code with Property 
(C,T>') of size q log l c l lo sl :D I. Instead, we show that it is possible to achieve a size of 
q K log l c I log l x> I for some positive constant k < 1. 

A straightforward generalization of Construction I to union of cosets of linear 
codes can be achieved as follows. Let C\ and T>\ be linear [n, fci,dci] and [to, 
such that j n G Ci and j m G Xfi. Let U C (F™/Ci) rep and V C (F™/£>i) rep . We 
consider the component codes C and T>\ where 


C = Ci + u, and V = |^J T)\ + v. 

u £U vGV 


Then the (to x n)-matrix code defined by 


(4.1) 


|J (J (Ci + u) ® (Vi + v). 

u £U vGV 


has Property (C',V). However, observe that the code has size 

\U\\V\q kl11 = g^Pi+log |W|+log |V| ^ 

while the sizes of C' and V are \U\q kl = q k i+ log l w l and |V|g Zl = g ; i+ log l v l respec¬ 
tively. Thus the size of the code obtained from (4.1) is less than g log l c l^sl 23 I = 

^Oi+log |W|)(Zi+log |V|) 
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4.1. Product Construction of Expurgated Codes. We improve the size 
given by (4.1) when the union of cosets of product codes has a certain structure. 
Specifically, we consider the instance where the cosets form an expurgated code. We 
describe this formally below. 

In addition to the codes 4, Pi, assume that 4 and V 2 are linear [n, 4, dc 2 ] and 
[m, 44> 2 ] codes such that 4 C C 2 and P 1 C P 2 . We consider nonlinear component 
codes that are obtained from expurgated codes C 2 \ 4 and P 2 \ P i- Our objective is 
therefore to construct an (m x n)-matrix code such that Property (4 \ 4, P 2 \ Pi) 
holds. 

Clearly, 4 \4 and P 2 \Pi are union of cosets of 4 and Pi with U = (4/4)re P \ 
{0„} and V = (P 2 /Pi)re P \ {Om} respectively. Then the construction described in 
(4.1) gives a code with size (g fc2_fcl — l)(q l2 ~ l1 — l)q kl11 ss q k ^~ k l+h-h+hih. qq ie 
distance of the product code is determined by the distance of the codes C 2 and P 2 . 

On the other hand, we improve this size via the following. 

Construction IA. Consider two intermediary codes 4 and P3 of dimensions 4 — 1 
and I 2 — 1 respectively such that 4 C 4 C 4 and Pi C P 3 c P2. Pick any 
u £ (4 \ 4) and v <E (P 2 \ P 3 ) and observe that 

C 3 + ucC 2 \Ci and P 3 + v C P 2 \ Pi. 

Applying Construction I to the cosets 4 + u and P 3 + v yields a matrix code 
(4 + u) ® (P 3 + v) with Property (4 + u, P 3 + v), and hence the Property (4 \ 
Ci, P 2 \ Pi). Furthermore, the size of this code is g^ 2-1 )^-!) anc [ j s significantly 
larger than the straightforward construction from (4.1). 

4.2. Binary Matrix Codes with Restricted Column and Row Weights. 

In this section, we apply Construction IA to obtain matrix codes with the additional 
property of bounded row and column weights. The motivation for studying such 
matrix codes arises from the application to coded modulation for power line commu¬ 
nication (PLC) channel. Consider a codeword N £ F™ xn of a matrix code. Each row 
of the matrix corresponds to transmission over a particular frequency slot, while each 
column of the matrix corresponds to a discrete time instance. Transmision occurs at 
the frequency and time slots corresponding to a one in the matrix. 

The different types of noises are as follows. Assuming a hard-decision threshold 
detector, the received signal (which may contain errors caused by noise) is demod¬ 
ulated to an output N £ F™ xn . The burst and random errors that arise from the 
different types of noises in the PLC channel (see [3, pp. 222-223]) have the following 
effects on the detector output. We denote the (z, j)-th entry of a matrix N by Nij. 

1. A narrowband noise introduces a tone at all time instances of the transmitted 
signals. If e £ [to] and e narrowband noise errors occur, then there is a set 
p e (M) 0 f e r0WS; suc p that Nij = 1 for i £ T, j £ [n\. 

2. Impulse noise results in the entire set of tones being received at a certain 
time instance. If e £ [ n\ and e impulse noise errors occur, then there is a set 
II £ (^) of e columns such that Nij = 1 for i £ [to], j £ II. 

3. A channel fade event erases a particular tone. If e £ [to], and e fades occur 
then there is a set T £ (^) of e rows such that Nij = 0 for all j £ [n]. 

4. Background noise flips the value of the bit at a particular tone and time 
instance. If e background noise occurs then there exists a set f l £ (["^[ m l) 
such that Nij = Nij + 1 , for all ( i,j ) £ Q. 
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We refer to [3] for an expanded description of the types of noise that are present in 
the power line channel. 

If any row of the codeword matrix N is an all-one vector then this row is not 
distinguishable from an all-one row introduced by the presence of narrowband noise. 
Similarly, an all-one column is not distinguishable from impulse noise. Additionally, 
the use of multi-tone frequency shift keying is adopted with the understanding that 
the energy is concentrated on a fraction of the available frequencies (see [7]). Thus, it 
is desired that every row and every column of the matrix contain at least a single one, 
but it should not be an all-one vector. This requires the use of codes whose codewords 
are matrices with restricted and bounded column and row weights. 

In particular, for the application to powerline communications we construct codes 
which are able to correct narrowband and impulse noise, random noise, and also 
simultaneously satisfy all of the following criteria (also see [7]): 

(Al) have positive rate, 

(A2) have positive relative distance, 

(A3) have efficient decoding algorithms, and 

(A4) have no restriction that the length of the code is at most the size of the 
alphabet. 

In the following text, we use the following lemma that was crucial in proving the 
so-called low symbol weight property (see [17, Proposition 1]) for g-ary affine codes. 

Lemma 4.1. Let C be binary linear [n,k,d] code such that (j ra ) C C. Then the 
codewords in C\ (j n ) have Hamming weight bounded between d and n — d. 

Proof. If the Hamming weight of any vector v is greater than n — d, then v + j„ 
has Hamming weight less than d. This is a contradiction. □ 

First, we illustrate via an example that the code obtained by straightforward 
expurgation does not satisfy the systematic property. 

Example 4.1. Let C = T> be the binary linear [4,3,2] code consisting of all even 
weight codewords. Observe that C \ (J 4 ) consists of six codewords of weight two and we 
are interested in constructing a (4 x 4 )-matrix code whose matrices have row weight 
two and column weight two. 

A naive approach is to look at the (3 x 3) information matrix and require all 
columns and rows to not belong to { 03 ^ 3 }. This approach fails as illustrated by the 
example codeword, 


(1 

0 

0 


0 

1 

0 

1 

0 

0 

1 

1 

V 1 

1 

1 

1 


which contains an all-one row even though each of the component codewords in the 
first three rows and columns have weight exactly two. 

On the other hand, consider the binary linear [4, 2, 2] code C 3 = {O 4 , J4, (1,0,1,0), 
(0,1,0,1)} and let u = (0,0,1,1). Then (C 3 + u) ® (C 3 + u) yields a (4 x 4 )-matrix 
code whose matrices have row weight two and column weight two. Furthermore, it is 
systematic of dimension four. 

On the other hand, it can be obtained via computer search that there are exactly 
90 matrices in C <S> C that have constant row weight two and constant column weight 
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two. An exhaustive computer search shows that there do not exist five coordinates 
where a subset of these 90 matrices is systematic. 

We proceed with the construction of matrix codes with restricted row and column 
weights. Let C, V be binary linear [n, k, dc], and [to, l, dp] codes respectively. Suppose 
(j„) C C, and (j m ) C D. Direct application of Construction IA yields a systematic 
binary (to x ro)-matrix code of dimension (k — 1 )(Z — 1 ) whose matrices have 

(i) row weight bounded between dc and n — dc, 

(ii) column weight bounded between dp and m — dp. 

In Example 4.1 we showed that this construction gives more desirable results 
and why naive methods of constructions do not work. Because of the narrowband 
and impulse noise present in the power line channel, we want codes with restricted 
column and row weights. The following proposition gives the condition under which 
the noises can be corrected. 

Proposition 4.2. Let C, V be binary linear [ n , k , dc ], and [to,/, dp] codes re¬ 
spectively. Suppose (j n ) C C, and (j m ) C T>. Then (C \ (j„)) <g) (T> \ (j m )) obtained 
using Construction I A yields a systematic binary (to x n)-matrix code of dimension 
( k — 1 )(Z — 1 ) whose matrices have 
(i) row weight bounded between dc and n — dc, 

(ii) column weight bounded between dp and m — dp. 

Furthermore, (C. \ (j„)) 0 (T>\ (j m )) is a subcode ofC®V, and hence, is able to correct 
ei\iBD narrowband errors and eiMP impulse noise errors, provided 

eiMP < dc, and cnbd < dp. 


Proof. Consider a code C' C C of dimension k — 1, and V CT> of dimension l — 1. 
Using Construction IA, we consider the cosets C' + u, and V + v where u £ C \ C and 
v £ V \ V. By Lemma 4.1, the weight of every vector in C + u is bounded between 
dc and n — dc- Similarly, condition (ii) holds. 

To correct cnbd narrowband noise errors and eiMP impulse noise errors, we use 
Algorithm 1 which first sets all narrowband noise and impulse noise errors to erasures, 
subsequently subtracts the coset leader, and then decodes the row and the column 
codes. In the absence of random errors, if the condition cimp < dc is satisfied, then 
the row code C' can correct all the corresponding erasures. Similarly, the column code 
T>' can correct all the erasures in each column if cmbd < dp. □ 

4.2.1. Optimality of the product construction. The affine codes obtained 
using the product construction Construction IA, are likely not optimal for large di¬ 
mensions of the matrix. Obtaining optimal codes which satisfy all the criteria (Al)- 
(A4) stated earlier in this subsection is still an open problem. Below, we show some 
examples of codes for which the construction is close to optimal. 

Example 4.2. Consider the first order Reed-Muller code with parameters C = 
[2 r , r + 1, 2 r-1 ]. The affine code (C \ (jn)) 0 (C \ (jn)) obtained by Construction IA is 
a (2 r x 2 r )-matrix code of dimension r 2 and Hamming distance 2 2r ~ 2 . This code can 
correct eiMP < 2 r_1 impulse noise errors and cnbd < 2 r ~ 1 narrowband noise errors. 


Before providing the next example, we recall a “Gabidulin Construction” from 
the thesis of the second author [13, Section 5.4] for square matrix codes. 
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Algorithm 1 : Decoder for Product of Affine Codes 

Input: detector output N € F™ xn , coset leader U 
Output: N'eC'o V' 


1 

2 

3 

4 

5 


/* Consider the narrowband noise as erasures 

for i G [to] do 

if N ij = 1 for all j £ [n] then 
] N i,j<—E for all j G [n] 

end 

end 


6 

7 

8 
9 

10 


/* Consider the impulse noise as erasures 
for j G [n] do 

if N ij G {l,e} for all i G [to] then 
] Nij <— s for all i G [to] 

end 

end 


11 

12 

13 

14 

15 


/* Subtract the coset leader from the nonerased coordinates 

for i G [m] , j G [n] do 
if N,;j ^ e then 
j ]N\; j ^— 1M 7; j — U 

end 

end 


16 Decode N to N' using a product code decoder 

17 return N' 


*/ 


*/ 


*/ 


Proposition 4.3 (Kiah [13]). Let d < n and let k = n — d+1. Then there exists 
a binary (2 n x 2n) matrix code of size 2 nk with constant column weight n that corrects 
eiMP impulse noise errors and cnbd narrowband noise errors provided that eiMP < n, 
cnbd < n, and 


eiMP I 

.“J + 


e NBD I 

. 2 J 


< d. 


The family of codes in the above proposition is obtained from Gabidulin codes [12], 
which are optimal rank metric codes, and can be explicitly written as follows. Let C 
denote a Gabidulin code, and let C* denote the matrix code obtained using Proposition 
4.3. Then we get, 


C* 


( M 

+ J 


M + J\ 
M ) 


M GC 


where J is the all-one matrix. We now proceed to provide an example similar to 
Example 4.2. 

Example 4.3. Consider the Gabidulin code with parameters [n = 2’ 1 ,1 ,d = 
2 r ~ 1 ]. Such a code has dimension nk = 2 r ~ 1 , and can correct the same number of 
narrowband and impulse noise errors that the code in Example 4-2 can correct. We get 
the following table comparing the dimensions of the two codes obtained from Example 
4-2 and Proposition 4-3 respectively. 
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r 

Dimension of codes from 
Example r 2 

Dimension of codes from 
Proposition 4-3, 2 r 1 

3 

9 

4 

4 

16 

8 

5 

25 

16 

6 

36 

32 

7 

49 

64 


Thus, beyond r = 6, the construction from the Gabidulin codes has better param¬ 
eters than the product construction. However, it must be noted that in the case of the 
Gabidulin construction, we have restricted the matrices in our matrix code to always 
be a square matrix. 

For fair comparison, we have taken the affine codes to also be formed of square 
matrices with the same number of rows and columns as in the Gabidulin construction. 
The product code in general does not have this restriction. This satisfies point (A4) in 
the criteria that we want to satisfy. Therefore, it is of interest (and an open problem) 
to obtain the true upper bound for codes which are optimal and which satisfy all the 
four criteria (A1)-(A4). 

5. Irregular Product of Affine Codes. The power line channel is known to 
be frequency selective (see [3]), i.e., the background noise in different frequency slots 
are of different intensities. Thus, it is of interest to provide constructions of codes 
that can provide different levels of error correction over different frequencies. Such 
codes can be constructed as product codes where the rows of the matrix correspond to 
different row codes. Such codes have been studied earlier as “generalized concatenated 
codes” or “multilevel concatenated codes” (see Blokh and Zyablov [ 6 ], Zinoviev [19], 
and Durner [9]). The row codes, which correspond to the row encoding, in these 
constructions are defined over an extension field of the field of the column code. 
As a result, although the resulting matrix code is linear over the smaller field, the 
rows do not in general belong to the row code. This makes it difficult to extend the 
construction to product of affine codes. 

Instead, we consider the case where the component codes for each row and col¬ 
umn are different. Although the application is only for binary component codes, 
we provide the general theory for q-ary component codes. Such product codes were 
termed irregular product codes and were studied by Alipour et al. [1]. Specifically, 
they demonstrated the following proposition. 

Proposition 5.1 (Alipour et al. [1]). Let Ci be a linear code of length n and 
dimension ki for i E [to] and T>j be a linear code of length n and dimension lj for 
j E [to]. Suppose that k\ < < • • • < k m and l\ < I 2 < • • • < l n . Then there exists a 

linear (in x n)-matrix code of dimension K, where 

n lj 

(5.1) I< < E E ma x{ki — j + 1,0}, where Iq = 0, 

1=1 i=lj- 1+1 

and every codeword N satisfies the properties that 
(i) the i-th row of N belongs to Ci for i E [to], and 
(ii) the j-th column o/N belongs to T>j for j E [n]. 

Furthermore, if Ci C C 2 C ••• C C m and T>i C T >2 C • • • C T> n , we achieve equality in 

(5.1) . 

The encoding algorithm of irregular product codes is described in [1], The en¬ 
coding procedure encodes the rows first and then encodes the columns. The encoding 
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assumes that the first fc.j coordinates of the *-th row can generate the remaining n — ki 
coordinates of that row, and that the first lj coordinates of the j-th column can gener¬ 
ate the remaining m — lj coordinates of that column. Since the generating coordinates 
of the code are present within the leading principal l n x k m submatrix, we have the 
following Lemma. 

Lemma 5.2. The leading principal l n x k m submatrix generates all the remaining 
coordinates of a codeword N in the irregular product code. 

We apply Construction I directly to Proposition 5.1 to obtain an irregular product 
of affine codes. 

Proposition 5.3. In addition to the conditions of Prop. 5.1, let j n G Ci for 
i G [m\ and j m G Dj for j G [n]. Let u= ( 0 km , a) G F” \ (J™ x Ci and v = (0; n , b) G 

*7\u; =1 2v 

Then there exists an affine (to x n)-matrix code of dimension K bounded by (5.1) 
and every codeword N in the code satisfies the properties that 
(i) the i-th row of N belongs to Ci + u for i G [m], and 
(ii) the j-th column of N belongs to Vj + v for j G [n]. 

If C\ C Cj C ■ • ■ C Cm and T>i C 2? 2 C • • • C T) n , we achieve equality in (5.1). 

For <7 = 2 , suppose there exist linear codes C and V such that (J”l 1 Ci C C and 
UHi 'C’i C V, respectively. For u G C \ (J)” x Cj and vGP \ Uj=i ^j> the weight of 
every row of any codeword is bounded between dc and n — dc, and of every column 
between dp and m — dp, where dc and dp are the minimum distances of C and V 
respectively. 

Proof. Let N be a codeword obtained by using the encoding described in [1]. We 
translate this codeword using the matrix 


Ofc m X/ n 

■T 

•l/„ a 

b T j k m 

bT j n-k m +j™-i n a 


where the vector a has length n— k m and b has length m—l n . We denote the codeword 
N by four submatrices, as 


N = 




where Ni is the l n x k rn leading principal submatrix that generates N 2 ,N 3 ,N 4 , by 
Lemma 5.2. The submatrix N 2 is of dimension l n x (n — k m ), N 3 is of dimension 
(to — l n ) x k m , and N 4 is of dimension (to — l n ) x (n — k m ). Denote the corresponding 
matrix from the coset code as 


N' = 




where Nj = N 2 + j^a, and N 3 = N 3 + b T jfc m . 

We need to ensure that the matrix N 4 obtained by encoding the rows of N ( 3 
by the row codes satisfies the condition that they belong to the row code, and also 

satisfies that they belong to the column codes in the submatrix 



. This can be 
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proved as follows. Let the generator matrices of Ci, i = l n + 1,..., m be given by 
the matrices Gj = \I ki |A'|Aj], where A\ has dimension fej x ( k m — ki) and A, has 
dimension ki x (n — k m ). Let b = (&; n+ i,..., b m ). The i-th row of Ng can be split 
into two parts, corresponding to the first two blocks of the generator matrix Gi as 
follows. For ! = l„ + l,...,mwe first write the i-th row of N 3 as Ng^ = (n^, n') where 
the first block has length ki and the second block has length k rn — ki. We obtain, 


N 3,i = ( n * + biiki, n- + bijk m -ki) = N 3ii + bijk m ■ 

Encoding the first block of this row Ng i with the generator matrix Gi gives the vector 

(w + bij ki )[I ki \A'i\Ai\ = (n, + 6 i j fei ,n i A' + bij km -ki, wAj + b,j n -k m ) 

= (n. i ,n-,N 4ji ) +b,j n 


where n, A' = n' is the second block of Na^, and n,Yl t = N 4 i is the *-th row of N 4 . 
The shift by the coset leader (Ofc m , a) results in the word 

(( n ii n 0 + frijfcm) N 4) j + bi] n -km + a) = (N' 3ii , N 4ji + bi] n -k m + a )- 

Thus, the matrix N 4 is given by the expression 

N 4 = N 4 + b T j„_ fcm + j™_ iri a. 


A similar argument shows that the submatrix 



satisfies the corresponding col¬ 


umn codes. 

The same argument as in the proof of Proposition 4.2 shows that the row and 
column weights are bounded when the conditions C\ C C 2 C • • • C T> rn C C, and 


2?i C V 2 C • • • C V m C V hold. □ 


6. Conclusion. We provide new constructions of systematic nonlinear product 
codes that are obtained by taking product of cosets of linear codes. The constructions 
have the property that every row and every column belongs to the row code and col¬ 
umn code, respectively. Subsequently, we show that it is possible to construct matrix 
codes with restricted column and row weights. Although the primary motivation for 
studying such matrix codes is for coded modulation over power line channel, the con¬ 
structions can potentially be adapted to address other problems where such codes are 
desired such as codes for memristor arrays and two-dimensional weight-constrained 
codes [15,16]. 
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